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This paper provides an overview of ongoing efforts to develop, evaluate, and validate 
different tools for improved aerodynamic modeling and systems analysis of Hybrid Wing 
Body (HWB) aircraft configurations. Results are being presented for the evaluation of 
different aerodynamic tools including panel methods, enhanced panel methods with viscous 
drag prediction, and computational fluid dynamics. Emphasis is placed on proper prediction 
of aerodynamic loads for structural sizing as well as viscous drag prediction to develop drag 
polars for HWB conceptual design optimization. Data from transonic wind tunnel tests at 
the Arnold Engineering Development Center’s 16-Foot Transonic Tunnel was used as a 
reference data set in order to evaluate the accuracy of the aerodynamic tools. Triangularized 
surface data and Vehicle Sketch Pad (VSP) models of an X-48B 2% scale wind tunnel model 
were used to generate input and model files for the different analysis tools. In support of 
ongoing HWB scaling studies within the NASA Environmentally Responsible Aviation 
(ERA) program, an improved finite element based structural analysis and weight estimation 
tool for HWB center bodies is currently under development. Aerodynamic results from these 
analyses are used to provide additional aerodynamic validation data. 
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I. Introduction 

H ybrid wing body (HWB) or Blended Wing Body (BWB) aircraft designs have been considered promising 
alternatives to conventional tube and wing aircraft designs due to their large potential fuel savings and 
increased aerodynamic efficiency/ Fuel burn reductions of 25% and higher have been published for some ultra-high 
capacity configurations carrying up to 800 passengers. Most notable is the work of Liebeck and his co-workers at 
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the Boeing Company. Their 450 passenger 
BWB-4501L design recently resulted in the 
8.5% scale X-48B flight demonstrator 
(Fig. 1)." 

A significant difficulty in dealing with 
HWB design optimization has always been 
the lack of a data base of known “flying” 
designs which may serve as calibration and 
validation points for optimization programs 
like FLOPS, especially when compared to 
the vast number of available tube and wing 
aircraft data. When transitioning from the 
conceptual to the preliminary design phase, 
the aircraft designer needs to be sure that the 
design chosen for further optimization is 
actually a viable design, and as a result of 
the lack of validation cases, the development 
of improved fidelity analysis tools becomes imperative already for the conceptual design loop. 

Therefore, in order to validate the projected fuel burn and noise reduction potential of HWB designs for NASA’s 
Environmentally Responsible Aviation Program (ERA), significant efforts have been put forward to develop 
advanced structural and aerodynamic analysis tools for HWB conceptual design optimization. Recently developed 
tools include finite element (EEM) based structural analyses to provide enhanced capabilities for HWB center body 
structural sizing and weight estimation. Aerodynamic methods improvement has been geared towards increased 
fidelity in-the-loop methods like enhanced panel codes and computational fluid dynamics (CED). 

While the use of CED methods to explore large design spaces is still limited by available computation power and 
reasonable turn-around times, several promising alternatives like enhanced panel methods with viscous drag 
prediction capabilities have become candidate tools for HWB conceptual design optimization. This paper describes 
recent efforts at NASA Langley’s Aeronautics Systems Analysis Branch towards the evaluation and validation of 
different aerodynamic tools ranging from panel codes to Euler and Navier-Stokes CED methods. Significant low- 
speed and transonic wind tunnel testing at the NASA Langley National Transonic Eacility (NTE) and the Arnold 
Engineering Design Center 16ft Transonic Tunnel (AEDC-16T) preceded flight testing of the X-48B and resulted in 
a wealth of available aerodynamic test data. Validation data for the present study is based on a June 2007 transonic 
test in the AEDC-16T on a 2% scale model of Boeing’s BWB-450-1L design. 

II. Aerodynamic Modeling and Analysis 

Eor aerodynamic modeling and validation, wind tunnel model data of the Boeing X-48B 2% scale model was 
used to generate geometry and input files for the different aerodynamic tools. Several different models were 
generated to export the X-48B geometry into required formats for the various analysis packages (Table 1). The 
greatest challenge was that each analysis tool requires different sets of input parameters, with different input file 
formats. 



Figure 1: Boeing’s X-48B Blended Wing Body flight 
demonstrator (Boeing picture). 


Table 1: Summary of investigated aerodynamics codes. 


Code name 

Code type 

Remarks 

Reference 

Vorlax 

Panel code, VLM 

Includes thick leading edge correction 

3 

Vorview 

Panel code, VLM 

Graphical front end for Vorlax 

4 

MSES 

2-D Euler 

Stripwise airfoil drag prediction 

5 

CPPAero 

Enhanced panel code, VLM 

Includes stripwise shock and viscous drag 

6 

Nastran 

Panel code, DLM 

Included in structural model for loads 

7 

CartSD 

CFD, 3-D Euler 

Used as scaling reference for full scale configuration 

8 

FunSD 

CFD, 3-D Navier-Stokes 

Used as scaling reference for full scale configuration 

9 
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A. Model Development 

As a first step, triangularized surface data 
of the X-48B wind tunnel model was 
imported into Matlab to generate suitable 
input files for the different analysis tools. 

Model data provided by the NASA Langley 
Configuration Aerodynamics Branch is a 
Tecplot data file of the 2% scale wind tunnel 
model of the Boeing BWB-450-1L with a 
wingspan of 4.80ft and a centerline chord 
length of 3.13ft. The model data is given by 
a set of unstructured surface grid data points 
from triangular panels (Fig. 2)}^ 

A Matlab routine was developed to 
generate hermite files for Vorview. To 
provide a cleaner surface for initial code 
validation of the vortex lattice panel 
methods, engine data was removed from the 

slices. The slice data was written to a text data file in the Vorview hermite format. To enhance portability of the 
model data to other codes available in the NASA Langley Aeronautics Systems Analysis Branch and to generate 
geometry files for the CFD codes, a VSP^^ model was developed from the extracted airfoil stack data. An overall 
modeling summary is shown in Fig. 3. 



Figure 2: Tecplot model of the X-48B 2% scale rotary wind 
tunnel model. 



Figure 3: Aerodynamic codes validation modeling overview. 
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B. HWB Nastran Model 

The HWB Nastran model has been developed in support of the HWB scaling study presented by Nickol^^ to 
investigate the viability of a hybrid wing body airplane for smaller commercial transport aircraft down to regional jet 
size (about 100 passengers). In an effort to improve previous FEM and FLOPS based analysis methods for weight 
estimation of the HWB center body, a full aeroelastic finite element model of an HWB primary structure has been 
implemented. To maintain the capability to interface with FLOPS, most of the structural parameters for finite 
element modeling of the HWB structure are based on the center body home plate analogy introduced by Nickol and 
McCullers.'^ 

As part of the aeroelastic validation of the structural model, Nastran aerodynamic loads have been compared to 
Vorlax and CPPAero results. The Nastran model is generated by exporting the VSP surface model and extracting a 
parametric representation of the HWB primary structure from the exported data, which is then used to build a 
Nastran bulk data file (Figs. 4 and 5). More details regarding the structural HWB finite element modeling in Nastran 
can be found in Ref. 14. 



Figure 4: Left: Exported VSP surface data with overlay of HWB primary structure. Right: Patran rendering 

of the generated Nastran bulk data file. 



Figure 5: HWB Primary structure with overlay of Nastran doublet lattice aerodynamic paneling. 
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III. Reference Wind Tunnel Data Set 

In June 2007, the X-48B 2% scale model of the BWB-450-1L was tested at the Arnold Engineering 
Development Center (AEDC) 16T wind tunnel. The test was called TE1038 under the AEDC designation system. 
As a result of this test, sufficient data was generated to develop a transonic aerodynamic database for vehicle lift, 
drag, and moment coefficients vs. angle of attack and Mach number. The data set used for this analysis was 
collected at Mach numbers ranging from 0.5 to 0.85 and MAC based Reynolds numbers of 10 million. Complete 
details of the AEDC T1038 test are given in Ref. 15. Results for vehicle C^, Cm, Cd, and drag polars vs. Mach 
number are shown in Eig. 6. 






Figure 6: Vehicle baseline aerodynamic coefficients from AEDC 16T transonic wind tunnel test. 


IV. Aerodynamic Validation Results: Panel Codes 

Vorlax is a generic vortex lattice panel code for subsonic and supersonic flow applications. To improve its user 
friendliness, the graphic front end Vorview has been wrapped around the original Vorlax solver. MSES was used for 
viscous drag predictions within Vorview. CPPAero is a VLM based aerodynamics analysis tool developed by Avid 
EEC in Yorktown, VA specifically for BWB configurations. The program uses Vorlax as a general VLM solver but 
has several additional, partly BWB specific features to account for compressibility and transonic effects, thickness 
effects, viscous drag, and BWB specific S&C issues. Through a built-in XML loader, an XML file can be imported 
into CPPAero in order to generate the input and geometry files required to run a CPPAero analysis. 
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A. Lift Coefficients 

Figure 7 shows lift coefficient comparisons between wind tunnel data, Vorlax, and CPPAero runs for the X-48B 
2% scale rotary wind tunnel model. Both low speed and transonic data sets agree exceptionally well with wind 
tunnel data in the linear region of the lift gradient. As a linear method, Vorlax does not predict the transonic lift rise 
nor flow separation at higher angles of attack. CPPAero ’s built-in corrections capture but overpredict the transonic 
lift rise up to pitch angles of stall onset. 





Figure 7: Lift coefficients vs. pitch angle and Mach number for the investigated panel codes. 

B. DragPolars 

Drag polars from all panel codes were compared to experimental data for the same Mach number range (Fig. 8). 
In general, panel codes only predict induced drag and portions of pressure drag related to airfoil camber. Therefore, 
different strategies have been used to improve drag predictions by including viscous drag calculations. The general 
approach for this process usually solves the linear VLM problem and then uses the VLM based pressure distribution 
to compute viscous drag contributions. 

CPPAero uses a zonal approach where the induced drag at a given subsonic/transonic condition is calculated 
using the incompressible drag polar shape as its basis. The wing Oswald efficiency factor is found from the 
incompressible Vorlax results. This is modified as necessary for the flow zones on each wing section (see Ref. 6). 

MSES uses 2D airfoil data to generate sectional viscous drag predictions which are summed up in the spanwise 
direction to provide the total viscous drag. Depending on the airfoil geometry and Mach/Reynolds number 
conditions, MSES has shown to have varying convergence characteristics. At lower Mach numbers, Vorview results 
with MSES viscous drag show excellent agreement with test data. For transonic Mach numbers, MSES convergence 
was sketchy, no reasonable result could be obtained for the Mach 0.85 condition. 
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Figure 8: Drag polars vs. Mach number for the investigated panel codes. 


C. Spanload Validation Using Nastran 

From the wind tunnel test data, only total lift, moment and drag coefficients were available. Spanload 
distributions for aerodynamics loads predictions were therefore validated by comparing vortex lattice Vorview 
results with spanwise lift 
coefficients obtained from 
Nastran doublet lattice aero- 
dynamics using the HWB 
structural model described in 
section II. B and Ref. 14. 

Figure 9 shows excellent 
agreement between both methods 
for several pitch angles ranging 
from 1 to 4 degrees. The Vorview 
spanload data is based on rigid 
lifting surfaces since Vorview 
does not employ an underlying 
structural model. Therefore, only 
rigid sectional lift coefficients 
from the Nastran solution are 

shown for the comparison. Figure 9: Vorlax vs. Nastran spanload comparison. 
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V. Aerodynamic Validation Results: CFD Codes 

Although not considered for direct use in a conceptual design and optimization loop due to computation power 
and turn-around time constraints, two CFD codes have been assessed in an effort to generate data for Reynolds 
number scaling from wind tunnel model dimensions to full scale HWB designs. CartSD is a high-fidelity in viscid 
analysis package for conceptual and preliminary aerodynamic design. It allows users to perform automated CFD 
analysis on complex geometries. The package includes utilities for geometry import, surface modeling and 
intersection, mesh generation, flow simulation and post-processing of results.^ 

Fun3D is an unstructured viscous CFD code solving the three-dimensional Navier-Stokes equations. The solver 
has been applied to a number of complex large-scale problems by groups at NASA, industry, and academia. 
Internally, the software has been used to study airframe noise, space transportation vehicles, flow control devices 
using synthetic jets, the design of wind tunnel and flight experiments, and so forth. Boeing, Lockheed, Cessna, New 
Piper, and others have used the tools for applications such as high-lift, cruise performance, and studies of 
revolutionary concepts.^ 

A. Lift Coefficients 

As expected, results from both CFD codes show good agreement with tunnel data for the lift coefficients. 
However, it was noted that CartSD significantly overpredicts the transonic lift increase at higher pitch angles. As a 
reference for the investigated panel codes lift coefficients from Vorlax are also included in the plots in Fig. 10. 






Figure 10: Lift coefficients vs. pitch angle and Mach number for the investigated CFD codes. 


8 

American Institute of Aeronautics and Astronautics 


B. DragPolars 

With Fun3D being the only viscous code in this assessment, the drag polar predicted by Fun3D shows the best 
agreement with test data. Pressure drag from Fun3D agrees well with Cart3D and the panel codes for all Mach 
numbers (Fig. 11). 






— AEDCT1038, Mach=0.85, Re=10M 
-•-Vorlax With Winglets, Mach=0.85 
-♦—FunBd, Mach =0.85, Rec=17.8M 
^^Fun3D, pressure drag only 
-A-Cart3D, Mach = 0.85 
— Cart3D, Mach = 0.85, with engines 


AEDCTL038, Mach=0.8, Re=10M 

-♦-Vorlax With Winglets, Mach=0.8 
Fun3d, Mach = 0.8, Rec=17.8M 
Fun3D, pressure dragonly 
-♦-Cart3D, Mach = 0.8 
^I^Cart3D, Mach = 0.8, with engines 


cD_tot 


cD_tot 


Figure 11: Drag polars vs. Mach number for the investigated CFD codes. 


VI. Drag Polar Discussion - CFD vs. Enhanced Panel Methods 

Subsonic drag predictions of both CFD and panel codes agree well with tunnel data. At transonic Mach numbers, 
the enhanced panel methods fail to predict viscous drag contributions accurately enough for aircraft performance 
optimization. Since both enhanced panel methods are based on stripwise viscous drag calculations, discrepancies in 
viscous drag appear to be mainly due to the difficulty in accounting for three-dimensional relief effects which 
become prominent at transonic Mach numbers. CPPAero accounts for 3-D relief based on the analytical 3-D relief 
relations introduced by Hoerner^^ and predicts total transonic drag trends more reasonably. As currently 
implemented, MSES viscous drag predictions are based on 2-D airfoil sectional drag with no 3-D relief. As a result, 
viscous drag is overpredicted at transonic Mach numbers. 

VII. Viscous Drag Scaling 

As shown in the previous sections, lift and moment coefficients from Vorlax agree very well with wind tunnel 
data and Vorlax spanload distributions have been verified by Nastran doublet lattice aerodynamics. In order to use 
Vorlax results for aerodynamic performance optimization, a more accurate method of viscous drag prediction was 
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derived with the goal of using induced drag numbers from Vorlax and adding viscous drag separately as a 
contribution to Cdo- This way, aerodynamic lift variations due to planform changes can be easily captured through 
Vorview’s geometry flexibility, while at the same time maintaining accurate viscous drag numbers without the 
necessity of excessive amounts of CFD runs. 

In general, viscous drag c^v can be derived from skin friction drag c^/by the well-known relationship 




wet yx 

/ ^ JJ HWB 

^ref 


( 1 ) 


In Eq. (1), jfnwB is a generic HWB form factor derived from CFD or wind tunnel data, while Sy^et and Sref are the 
wetted and reference areas of the HWB, respectively. 

A fundamental relationship for the skin friction drag C/ vs. Reynolds number has been presented by Hoerner as^^ 

c^=i^/Ref'" (2) 


Ref. 16 indicates that for a Reynolds number range between lO’ and 10®, the values m = l and K = 0.030 may be 
used. The resulting relationship between Reynolds number and skin friction coefficient C/is shown in Fig. 12 as a 
blue line. To account for compressibility effects, additional corrections to Hoerner ’s skin friction drag relations have 
been applied based on free stream Mach number^^'^^, where 


c 


fcom 



(3) 


with 


r = l + 0.155M^ 


(4) 


To check the scalability of the presented approach for HWB viscous drag predictions, additional CFD runs were 
performed at different Reynolds numbers from subscale Reynolds numbers to full scale cruise flight conditions for 
the Boeing BWB450. CFD results match the Hoerner equation extremely well across the entire Reynolds number 
range from 5M to 200M, i.e. covering both the X-48B 2% scale wind tunnel model and BWB450 full scale cruise 
conditions (see dashed lines in Fig. 12). Additional available CFD based viscous drag data for the BWB710 follows 
the same trend and is also labeled in Fig. 12. 
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Figure 12: Skin friction drag and viscous drag coefficients vs. Reynolds number from CFD data and wind 

tunnel testing. 
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Drag predictions from the 2% scale AEDC T1038 wind tunnel test show that viscous drag is slightly 
underpredicted by Fun3d. This is a common issue in viscous (and therefore total) drag prediction and may be 
attributed to wind tunnel installation effects, CFD turbulence models, differences in wind tunnel model vs. CFD 
model configuration etc. 

HWB form factors for Eq. (1) have been derived from both CFD viscous drag results and wind tunnel test 
data. Total drag coefficients vs. pitch angle as well as drag polars for both form factors are shown in Figs. 13 and 14. 






alpha (deg) alpha (deg) 


Figure 13: Total drag coefficients vs. pitch angle for Vorlax with viscous drag (form factors derived from 

CFD and wind tunnel data). 


As expected, total drag predictions using CFD based form factors slightly underpredict wind tunnel drag 
measurements. Vorlax pressure drag with wind tunnel based viscous drag predictions shows excellent agreement 
with tunnel data even at transonic Mach numbers. While total drag predictions depart from tunnel data at higher 
pitch angles, predictions match tunnel data extremely well in the pitch angle and Mach number range relevant for 
aerodynamic performance optimization at cruise, as well as low speed take-off and landing conditions. 

The Hoerner approach for predicting viscous drag for HWB configurations works well because the form factors 
derived from wind tunnel or CFD data are rather insensitive to planform geometry variations within this aircraft 
category. This method of calculating viscous and total HWB drag for aerodynamic performance optimization has 
been successfully applied by Nickol in his HWB scaling study. 
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Figure 14: Drag polars for Vorlax with viscous drag (form factors derived from CFD and wind tunnel data). 


VIII. Conclusions 

Lift and moment predictions from both CFD and panel codes are suitable for loads calculations and HWB sizing. 
As expected, at transonic Mach numbers, fully viscous Navier-Stokes based CFD tools like Fun3D outperform 
enhanced panel methods and Euler codes in accurately predicting total drag polars for performance optimization, 
mainly due to the difficulty of accurately accounting for three-dimensional relief effects or viscous boundary layer 
effects. As of now, all of the investigated enhanced panel codes failed to provide accurate viscous drag predictions. 
Results obtained from CPPAero demonstrate that this problem may be overcome by implementing analytical 3-D 
relief procedures. However, while trends are predicted correctly, increased accuracy will still require tuning of the 
respective program parameters with wind tunnel or CFD data. 

The use of a Reynolds number and Mach number based scaling scheme derived from the Hoerner relationships 
for predicting skin friction and viscous drag so far has shown the best results for viscous drag predictions apart from 
using fully viscous CFD analyses. The Hoerner approach works well for HWB configurations since form factors for 
viscous drag calculations are deemed to be rather insensitive to HWB planform geometry variations due to their 
overall aerodynamic simplicity and effectiveness. After all, for this very reason, HWB airplanes are being 
considered as alternatives to conventional tube and wing transport airplanes.^ 
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